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Abstract 

We consider the problem of determining a pair of functions (u, /) satisfying 

^ I the heat equation ut — Au = ip{t)f{x, y), where (x, y) € = (0, 1) x (0, 1) and the 

■^ ' function if is given. The problem is ill-posed. Under a slight condition on ip, we 

,J^ . show that the solution is determined uniquely from some boundary data and the 

C^ I initial temperature. Using the interpolation method and the truncated Fourier 

series, we construct a regularized solution of the source term / from non-smooth 

data. The error estimate and numerical experiments are given. 
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1. Introduction 



O I Let T > and let Q = (0,1) x (0,1) be a heat conduction body. We consider the 

problem of determining a pair of functions {u, f) satisfying the system 



' Ut- Au = ip{t)f{x,y), 

u^{0, y, t) = u^{l, y, t) = Uy{x, 0, t) = Uy{x, 1, t) = 0, 

uil,y,t) =0, 
^u{x,y,0) = g{x,y), 

for {x,y) E fi, t G (0)^)) where g G L^{Q) and (f G L^(0,T) are given. 

This is a case of the problem of finding the source F{^,t) satisfying the heat equation 

Ut — Au = F, 

where ^ is the spacial variable. This inverse source problem is ill-posed. Indeed, a 
solution corresponding to the given data may not exist, and even if the solution exists 
(uniquely) then it may not depend continuously on the data. Because the problem is 
severely ill-posed and difficult, many presumptions on the form of the heat source are 
required. Roughly speaking, the function F can be approximated by 

N 

F(e,t)^5^V^„(t)/„(0. 

n=0 



For simplicity, one reduces this approximation to its first term 

where one of the two functions Lp and / is given. Many authors considered the uniqueness 
and stabihty conditions of the determination of the heat source under this separate form 
[21 [HI [HI [15], [16], [7l [HI H] . In spire of the uniqueness and stabihty resuhs, the regularization 
problem for unstable cases is still difficult. For a long time, it has been investigated for 
the heat souce which is time-depending only [I2l[9l[5] or space-depending only p!l[T3l[5]. 
Recently, the regularization problem for the heat source F{^,t) = ip{t)f{^), where ip 
is a given function, was regarded for one-dimensional case [10] and two-dimensional 
case [11]. However, these authors needed in addition an essential datum, that is the 
the final condition u{^,T). Although this condition is unnatural, it gives an explicit 
representation of the solution as the inverse Fourier transform of a known term, and 
hence ones could use the truncated integral method to construct a regularized solution. 

In the present paper, we consider a similar problem to [TOl[Tl], but the final condition 
is removed completely. Moreover, the overspecified condition, i.e. u{l,y,t), is slighter 
than this one in flT] and is almost optimal to still hold the uniqueness of the solution (see 
Remark 3). To our knowledge, no explicit form of the solution of system ([T|) is available, 
and hence it is not easy to solve the problem although one has exact data. Of course, 
the problem with approximate data is even more difficult because of the ill-posedness. 

Under a slight condition on (/?, we shall use the variational method and some proper- 
ties of analytic functions to show the uniqueness of the solution. In particular, this result 
makes a regularization theorem of pllj trivial (see Remark 2). In spite of the uniqueness 
result, the problem is still ill-posed. We mention that the existence problem of a solution 
is not considered here. Instead, we shall assume that there is a (unique) exact solution 
corresponding to the exact data, and our aim is of constructing a regularized solution 
from approximate data. Using the interpolation method, we shall seek the coefficients 
of the Fourier series expansion of the source term / and then construct a regularized 
solution by the truncated Fourier series. The error estimate between the regularized 
solution and the exact solution is of order (ln(£:^^))~^, where e is the error between the 
given data and the exact data. We also note that we shall concentrate only on finding 
the source term / because we shall get a classical heat problem as soon as we know this 
function. 

The remainder of the paper is divided into four sections. We shall introduce some 
notations and state main results in Section 2. After that, we shall prove the uniqueness 
result in Section 3 and the regularization result in Section 4. In Section 5, we shall show 
how our method can be numerically implemented and give two examples to illuminate 
its effect. 

2. Notations and main results 

By variational method, we have the following formula to reconstruct the solution of 
the system (1). 

Lemma 1. Assume that {u J) e {C\[0,T]; L\n)) n L^{0,T; H^{n)), L^{n)) is a solu- 



tion of the system (1). Then for all {a, n) & C x Z we get 

g-{Q! -n n )T / iii^x, y ,T) cosh.{ax) cos{mry)dxdy 

— / g{x, y) cosh(ax) cos{mry)dxdy 



n 

T 

g-(a -n vr )*(^(^^^(^^^ / f (^x, y) cosh{ax) COS {mry)dxdy . 
n 

From Lemma 1, we introduce some useful notations. For w G L^{Q), (p G L^(0,T) 
and a, l3 E C, put 

G{w){a,(3) = w{x,y)cosh{ax)cos{Py)dxdy, 

Q 
T 



D{^){a,P) = / e-(" -/^ )V(0^^ 







^(-)(«'^) ifDM(«,/^)^0, 



0, ifZ}(v9)(a,/?) = 0. 

Note that for w G L^(r2) and integers m, n, 

G{w){irmi,mx) = w{x,y)cos{nnrx)cos{mry)dxdy. (2) 



n 



Because < \/K(m, n) cos{rmTx) cos{miy) > is an orthonormal basis on L^(i7), where 

I J m>0.n>0 



1, if m = n = 0, 
K(m, n) = ■^ 2, if m > n = or n > m = 0, 
4, if m > and n > 0, 

we have the following representation 

w{x,y) = y^ K{m,n)G{w){imn,mr)cos{mnx)cos{n7ry). (3) 

m,n>0 

This formula allows us to recover / from G{f). From Lemma 1, if (a^ — n^n^) > 
is large and \D{ip){a,n7T)\ is not so small then (?(/)(«, nyr) can be approximated by 
H{Lp,g){a,mv). To control |D(v9)|, we need the following condition (H) on (f. 

(H) There exist To G (0,T], ^ > and A > and such that either Lp{t) > At^ for a.e 
t G (0,To), or ip{t) < -At^ for a.e t G (0,To). 



Remark 1. The class of functions satisfying (H) is very broad. This condition holds 
with respect to 9 = 0, for example, if ^p is continuous att = and ip{0) 7^ 0. To compare, 
we refer to the condition tp G C^[0, T] and ip{0) 7^ in [I^, 



Under the condition (H), we will obtain the uniqueness of the problem ([T]). 

Theorem 1. Assume that g G L^{Q) and (f satifies (H). Then the system (Ql) has at 
most one solution (uj) m \c^{^,T\,L^{^))^L^{{S,T;E\^)),L^{^)) . 

In spite of the uniqueness, the problem is still ill-posed and hence a regularization is 
necessary. Generally, the main ideas of the regularization are divided into three steps. 
For each integer rz, we first approximate G{f){a,mi) by H{Lp,g){a,nn) for some real 
numbers a. In the next step we recover G{f){z,mr) when z is in a ball of the complex 
plane. Finally, we use a truncated series from the formula ([3]) to construct the regularized 
solution. 

For each integer n, we shall use the Lagrange interpolation polynomial to handle the 
key point of recovering G{f){., rni). Recall that if ^4 = {xi, ..., Xp} be a set of jo mutually 
distinct complex numbers and whe a, complex function then the Lagrange interpolation 
polynomial L[A] w\ is 

V / 



L[A-M{z) = Y.[\{ir^]w{x,) 



Xj Xk 



Now we are ready to state the regularization result. 

Theorem 2. Let go G L^{Q) and let ipo G L^{0,T) satisfy (H). Assume that {uq, fo) G 
{C\[0,T];L\n))nL\0,T;H^{n)),L^{n)) IS the exact solution of the system (1) cor- 
responding to the exact data go and ipo. Let e > 0, ipe & L^{0,T) and Qs G L^(0, 1)) such 
that 

The regularized solution f^ is constructed from (/^^ and g^ as follows 



r^ e zn 



In(e-i) In(e-i) 



50 50 

B{r,) = {±(4r, + j)|j = l,2,...,20rJ, 
Feim,n) = L[B{r^);H{ip^,g^){.,n7c)]{im'K), 
fe{x,y) = y^ K,(m,n)Fe(m,n) cos{m7Tx) cos{n7Ty). 

0<Tn,n<rir 

Then /, G C^(R'^) and lim f, = /o m L'^iVt). Moreover, if fo G H\n) then lim f, = 

e^0+ e^0+ 

/o in H^{Q) and there exists a constant £0 > depending only on the exact data such 

that 

50 

ll/e - /o|lL2(f7) < ^1^(^-1) WfoWmin) 
for all e G (0,£o)- 



3. Uniqueness 

Proof of Lemma 1. 

Proof. Getting tlie inner product in L^(f2) of the first equation of the system (1) and 
W{x, y) = cosh(ax) cos{n7vy), then using the integral by part we have 



d_ 
It 



uWdxdy - (a^ - n^vr^) / uWdxdy = (p{t) / fWdxdy. 



d_ 
It 



Next, we multiply the latter equality with e *^° "■ '^ )* to get 

g-(a2_„2,2), f^,^^^y\ _ g-(a2-n2.2)i^^^^ f fWdxdy. (4) 

n / n 

Finally, integrating (jl]) with respect to t from to T we obtain the desired result. D 

Now we consider some properties of the functions G{w) and D{ip). 
Lemma 2. Let w G L^{Q) and n & Z. Then G{w){.,mr) is an entire function and 

\G{w){z,mT)\ < e'"'' ||t«||ii(Q) , ^z e C. 

Moreover, if w & L'^{Q) and w ^ then there exists an integer n such that 

limsup — ' ^ ^^ ' > -1. 



1 ^ + 00 



Proof. For each integer n, put $n(-z) = G{w){z,n7i). Note that z ^— i> $n(^-2) is the 
cos- Fourier transform of the function 



1 
x I— > / w{x,y) cos{mTy)dy, x G [0,1]. 





and hence ^niiz) as well as $n(^) are entire functions. Moreover, 



I'^nl-z)! < / \w{x,y)cosh.{zx)cos{mTy)\dxdy 
n 



Now assume that w ^ 0. Note that 



d 
— $„(im7r) = / ixw{x,y)sin{mnx) cos{nny)dxdy. 



a 



and {sin(m7ra;) cos(?T-7r?/)}m>i^„>o is a orthogonal basis on L^(f2). Therefore, $„ is not 
constant for some integer n, and hence M<^^{r) > 1 for r > large enough, where 
Mq,^{r) = max|$n(z)|. We shall show that 

\z\=r 



r ln|'^n(r)| ^ 
hmsup > — 1. 

r— >+oo 1" 



Of course, it is sufficient to consider the case linisupln |$„(r)| < 0. For r > large 



r— >+oo 



enough, since In |$(r)| < and 1 < M^^{r) < e^ W'^j^Wl'^iq)^ we have 



ln|<l>„,(r) 



> 



ln|$nW 



r + ln||w7||^i(f^) In M$„(r) 



Moreover, according to Beurling theorem (see, e.g., ^, Section 6.1, page 40) we get 

> -1. 



,. In $„ r 

lim sup - — — — —^ 
r^+oo In M$fr) 



D 



This implies the desired result. 
Lemma 3. Let tp G L^{0,T) and {a,n) & Rx Z. Then 

\D{^){a,mr)\ < \\^\\l^o,t) 
when (a^ — n^yr^) > 0. Moreover, if ip satisfies (H) then 

hminf ia'-n\'f^'\D{^){a,nrr)\>0, 

where 9 is as in [H) corresponding to ip. 

Proof. The first inequality is obvious. Now assume that ip satisfies the condition {H) 
corresponding to Tq, 9 and A. We shall prove that 



liminfA^'+i 

A— ►+00 



e'^^ip{t)dt 



>0, 



which will imply the desired result by choosing X = a — n % . We have 



A' 



'+1 



e'^^ip{t)dt 



> A^+i 



To 



e~^'ip{t)dt 



A' 



'+1 



To 



To 



> AA^+1 / e-^H'dt - A^+ie^^^° 



e'^^ip{t)dt 



\'^\\l^o,t) 



for all A > 0. Since lim (A^+^e ^^°) = 0, it is sufficient to show that 

liminf ^^(A) > 0, 

A— >+oo 



(5) 



where 

To 



me{X) = \'+^ I e-^H'dt, VA>0. 



Using the integral by part we get 

^>e+i{\) = -{XToY^'e-'^o + {9 + 1)M/,(A). 

Therefore, it is enough to prove ([5]) for all 9 G [0, 1). Indeed, by direct calculus we obtain 
lim \l'i(A) = 1, lim \I/2(A) = 2, and by Holder inequality we get 

A— >+oo A— >+oo 

(^e(A))^ (*2(A))^ > *i(A), \/9 e [0, 1). 
Thus liminf ^6»(A) > 2^^^ for all 9 G [0, 1), and the proof is completed. D 

A— >+oo 

Proof of Theorem 1 

Proof. Let (wi, /i) and (m2, /2) be two solutions of the system (1). Put u = ui — U2 and 
/ = /i — /2. We have to show that (u, /) = (0, 0). 

Assume that / 7^ 0. It follows from Lemma [1] that, for all (a,n) & Rx Z, 



-(a^— n^TT^ 



^^G{u{., ., T)){a, uTv) = D{if){a, mi).G{f){a, mr). (6) 



Due to Lemma [21 there exist an integer n and a sequence of positive numbers {am}, 
which depends on n, such that lim a^ = +00 and, for all m > 1, 

m—*oo 

|G(/)(a™,n7r)|>e-2"'". 
Using Lemma [2] again we get 

|GK.,.,T))(a™,n^)|<e"-||n(.,.,T)||^,(^). 

Moreover, according to Lemma [3l there exist constants 9 > and Co > 0, which depend 
on Lp, such that 



\D{(p){am,mr)\ > 



Co 



[al - nH^y+^ 



for all m large enough. Therefore, from ([6]) we obtain that 



for all m large enough. This is a contradiction. 

Thus / = 0. Hence, the equality (jl]) in the proof of Lemma 1 reduces to 



d 



It follows from «(.,., 0) =0 that u = Q. D 



Remark 2. In Theorem 2 of ITl^ . the authors imposed that u{.,.,0) = and ip < 
—Co a.e. for some positive constant Cq. However, due to our uniqueness result, these 
requirements are so strict that the problem has only the trivial solution. 

Remark 3. In comparison to Ol]/ . here the datum u{x,l,t) is omitted. Moreover, if 
u{l, y, t) is omitted too then the uniqueness of the solution may not hold even if the final 
temperature is given. For example, the system 

Ut — /S.U = (TTCos^nt) + 27r^ sin(7rt))/(a;, y), 

M^(0, y, t) = 'Ua;(l, y, t) = Uy{x, 0, t) = Uy{X, 1, t) = 0, 

uix,y,0) = 0,u{x,y,l) = 0, 

has the trivial solution [u, /) = (0, 0) and (at least) a non-trivial solution 

u{x, y, t) = sin(7rt) cos(7rx) cos(7r|/), f{x, y) = cos(7rx) cos^Tvy). 

Thus the overspecified condition u{l,y,t) is essential. 

4. Regularization 

To construct the regularized solution, we first want to solve the problem of recov- 
ering the entire function G{f){.,mT) for each integer n. The key tool is the Lagrange 
interpolation polynomial. 

Lemma 4. Let r > 1 be an integer and B^ = {±(4r + j)|j = 1, 2, ..., 20r}. Let w and 
w be two even complex function such that w is an entire function and \w{z)\ < Ae'^' for 
all z & C , where A is independent on z. Then 

sup \w{z) - L{Br; w){z)\ < Ae~'' + 20re^^'^ sup \w{z) - w{z)\. 

|2|<7rT- zeBr 

Proof. Fix z E C, \z\ < nr and denote Zj = ir + j for each j = 1, 2, ..., 20r. We shall use 
the triangle inequality 

\w{z) - L{Br] w){z)\ < \w{z) - L{B,- w){z)\ + |L(5^; w - w){z)\. 

We first estimate \w{z) — L^B,,.; w){z)\. Let 7 = {z G C, \z\ = 45r}. Using the residue 
theorem, we get the Hermite's remainder formula 

-. „ /^x 20r 2 2 

Therefore, 

{, /pN I 20r I 2 _ 21 "j 



For (^ G 7 we have |tf(OI ^ Ae'^^^, \^ — z\ > (45 — 7i)r and 

We shall show that 

n ("^)' + ^^ .45-vr 

11 (45^)2 _^P 45 ^ ' ^^ ^'^'- ^y^ 

We can check ([9]) by direct computations for r = 1,2, ...,54. Now we consider when 
r > 55. Since the real function 



is increasing and concave in [0, (24r)^], we can apply Jensen's inequality to get 

20r 5 / 4A,T \ ^ / 1 ^'^ \ 

E''(s') = E E ''(-DU^-'E'' i E -n 

J=l fc=l yj=4(fc-l)r+l J fc=l y j=4(fc-l)r+l / 

= 4r^^( M6P + 16A; + — Jr2 + (4A; + 2)r + -J 

< 4rV^( (l6A;2 + 16A; + — ) r2 + (4fc + 2)— + 
fc=i VV 3 7 55 

= 4r^^('('l6A;2 + — fc + 

k=l ^ ^ 

5 

= 4r Y^ In 



6 X 55^ 
884, 32487' 



fc=i 



6050 

vr^+(l6P + tfc+lif) 
45^-(l6P + t^+lif) 



/45-7r\ 
= 4r X (-11.51809713) < -46r + In I , r > 55. 

Thus ^ holds for all r = 1, 2, .... From ([8]) and ©, we reduce © to 

\w{z)-L{Br;w){z)\ <Ae-\ Vr = l,2, ... (10) 

We shall now estimate 1-^(1?^; w — w){z)\. Since 

20r / ^2_^2\ 

L{Br] w~w){z) = ^\W 2_ 2 (^(^i) ~ ^(^i)) ' 
we obtain 

I 1 9 9 

\L{Br]w-w){z)\<2Qra sup TTft^^' (l^) 

l<i<20r^_^^. \Zj - Z^\ 



where a = sup \w{zj) — w{zj)\. We have 

l<j<20r 

Mi 1^^' - ''\ l}j ^'^ + ^^)^^- 1^^- " ^^1 ~ ti 1^^- " ^^1 

20r 

^2^" (4r + 2).(4r + 3)...(24r) 



< 



U\zj-Zk\ (j-l)!(20r-j)! 



^ (4r + 2).(4r + 3)...(24r) _ 

(10r-l)!(10r)! -^ "^IH- 

Note that J(l) < e^^ and 

J{r + 1) (24r + l).(24r + 2)...(24r + 24) 



Jir) (4r + 2)...(4r + 5) x (lOr) x [(lOr + l)...(10r + 9)] x (lOr + 10) 

24^^ 25 

< < e 

44 X 1020 

Here, the latter inequahty can checked by direct expansion 

24^^ X (4r + 2)...(4r + 5) x (lOr) x [(lOr + l)...(10r + 9)]^ x (lOr + 10) 
-4^ X 10^° X (24r + l)(24r + 2)...(24r + 24) 

to get a polynomial of degree 23 in term of r, which is obviously positive for r > 1. Thus 
J{r) < e^^^ for each integer r > 1. Therefore, we can reduce fllip to 

\L{Br;w- w){z)\ < 20re^^'a. (12) 

From (TTO!) and (TT2|) we have the desired result. D 

In our application of Lemma 4, we shall choose w = G'(/o)(., nyr), an unknown 
function, and w = H{ip^, gi;){.,mi'), which is known from data. We have the following 
error estimate between these functions. 

Lemma 5. Let uq, fQ,LpQ,gQ,Lp^,g^,r^,B{r^) he as in Theorem 2 and 6 he as in the 
condition (H) corresponding to ipo. Then there exists £1 > depending only on the 
exact data such that 



\Gifo)ia,nn)-Hiip,,g,)ia,nn)\ < {He-')) 



u^^e+5 I 



e' '£. 



provided e € (0, £1), < n < r^ and a G B{rs). 

Proof. Note that if a G B{r^) then 4r^ < \a\ < 24r^. Hence, for 5 > small enough one 
has 

^^1!^ < a^ ^ n^n^ < {\n{e-'))\ 
Thus, according to Lemma [3l there exists C{ipo) > depending only on (fo such that 

10 



and consequently, 



\D{ipe){a,mr)\ > \D{ipo){(^,n7r)\ - \D{ipo){a,nn) - D{ips){a,mr)\ 

_UN -2(0+1) 



It follows from Lemma [1] and Lemma [21 for e > small enough, that 



^ > 2 H^'')) 



|G(/o)(a,n^)-i/(^o,<7o)(a,r27r)|=e-("'-"''^')^ 
el"l|l«o(.,.,T)||^. 



G{uo{.,.,T)){a,mT) 



D{ipo){a,n7T) 



< g-(a2-„2^2)r 



(^) < le.el-l (ln(5-i))''^+' 



• (ln(e-i))-^(^+^) 
Moreover, for 5 > small enough, 

\H{ipo,go){a,n-K) - H{ipe, ge)ia,mr)\ 
G{9oi-,-,T)) (0,7171) G{gei.,.,T)){a,n7r) 



< 



D{ipo){a,n7T) D{ip,){a,mT) 

\G{go{.,.,T)){a,n-K)\ . \D{ip^){a,n7r) - D{ipo){a,nTc)\ 
\D{(po) (a, mr)D{ip^) (a, nn) \ 
\D{ipo){a,mr)\ . \G{gei., .,T)){a,mT) - G{goi., .,T)){a,mT)\ 



+ 



\D{ipo){a, mr)D{ip^){a, mr)\ 



< 



cl"^! II /I II c- 



(ln(e-i))-'('+') .2 (ln(£-i))- 



■2(6»+l) 



+ 



IV'o|Ili(o,T) • 



e'^le 



(ln(5-i))~'^'^"*''^ .2 (ln(£-i)) 



nN-2(9+i) 



^ 1 /l / -l\\49+5 Ul 



D 



The desired result follows the two latter inequalities and the triangle inequality. 
For w G L'^{Q) and an integer M > 1, put 

^Miw){x,y) = 2, fi(rn,n)G{irmT,mT) cos{rmTx) cos{mTy). 

0<m,n<M 

The formula ([3]) means lim TmIw) = w in L'^{Q). However, to prove Theorem 2 we 

Af^+oo 

shall need a sharper estimate for the remainder of the Fourier series. 
Lemma 6. If w & H^{Q) then lim TMiw) = w in H^{Q) and 

A'/^+oo 



1 



Mhu 



IIFmH - ^|lL2(f^) < ;^p;^— yy ii-'ii//i(q) 
Proof. From [2], using the integral by part we get 

7tmG{w){irmr,mT) = — / Wx{x,y)sm{m-iTx)cos{n-iTy)dxdy, 



TxnG{w) {imn, nn) 



Wy{x, y) cos{rmix) sm{mTy)dxdy, 



11 



and hence 

7r^(m^ + n'^)n{m, n) \G{w){im'K , mi)\ 

K{m,n) I / Wx{x,y)sm{rmTx)cos{nny)dxdy 

K{m, n) I / Wy{x, y) cos(m7rx) sm{mxy)dxdy 
\n / 

Since {sin(m7rx) cosiwKy)} ^^-^ ^>g and {cos(m7rx) mi(n-ny)} ^^^ ^^^ are orthogonal bases 
on L^(r2), it follows from Parseval equality that 



^ K{m^n) I / Wx{x,y)s\.\i{mTxx)cos{mTy) 



I 11^ 



m,n>0 



^ K{m,n) I / tfy(x, y) cos(m7rx) sin(n7r?/) 



I 11^ 

m,n>0 



Thus 



y 7r'^{ni^ + n'^)K{m.,n)\G{w){irn7r,mr)\ 

m,n>0 

— II l|2 I II ||2 ^ II ||2 

Using Parseval equality again, from the latter inequality we have 

\\^m{w)-w\\12^^) 
= 2, K^mjU) \G{w){i'miT,mT)\ 

max{m,n}>M 

^ 7r2(M + l)2 5Z ^^(m2 + n^)/t(m,n)|G'H(2m7r,ri7r)|^ 

max{m,n}>M 
1 2 



and 



irMH-wllii(f^) 



= 2, (l + 7r^(m^ + n^)) /«(m, 77.) |G'(w)(2m7r, ?7,7r)| ^0 

raax{m,n}>M 

as M ^ +00. D 

Proof of Theorem 2. 

12 



Proof. We shall first get the error between Tr^ifo) and fs, and then use the approxima- 
tion between r.r^(/o) and /o to get the results. Note that 

„.^,,. , F,{m,n) = L[B,;H{ip„g,){.,nn)]{irmr),ifO<m,n<r^, 

G{fe)[irmi,mT) = { . 

U, otherwise. 



Let us consider < ?72,n < r^. By Lemma [2l G{fQ){.,nn) is an entire function and 



|G(/o)(-2, n7r)| < ll/olliim) e'^', z G C. Moreover, for £ > small enough, Lemma [5] gives 



sup \G{fo){a,nn) - H{cp,, g,){a,n7T)\ < ilnie-')y'^'e'''^e. 

a£B{rs) 

Therefore, it follows from Lemma H] that 

\G{fo){im7r,mv) - G{fe){irmi,mi)\ 
= \G{fo){im7r,nn) - L[Be] H{ipe, ge)i-,mT)] (imn)] 

< ||/o|lii(f,)e-^^ +20r,e25^'^ {He''))''''' e''^^ e , 0<m,n<r,, 

for e > small enough. Note that e^'^ < ee~'^^^, the above inequality reduces to 

\G{fo){innT, rnv) — G'(/e)(«m7r, rnr)] < e^'^', < m, n < r^, 
for e > small enough. Therefore, 

W^rAfo) ~ fe\\m(n) 
= y, [l + 7i'^{rn'^ + n^)) K{'m,n)\G{fo){irnTT,nTT) — G{fe){irn7T,nn)\ 

0<m,n<r^ 

<(l + ^j2(l + 27rV,2) 45^/51 <eV26 

for e > small enough. 

Thus lim (r^_^(/o) — Z^) = in H'{Q). Hence it follows from lim Tr^lfo) = /o in 

L\n) that"" lim f, = /o in L^(n). 

Now assume in addition that /o G H'(Q). Then Lemma [6] leads to lim TrJfo) = /o 

in H'(Q) , and hence lim /^ = /o in H'(Q). Moreover, using Lemma 6 again we get 

Wfe - fo\\L2(n) < W^rAfo) - fe\\L2(n) + W^r.ifo) - fo\\L2(n) 

< pV52 , ^ II f II 

^ ^ ^— ^ ll/o|l/fi(C) 



< 



7r(r^ + 1) 



for £ > small enough. D 

5. Numerical experiments 
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In this section, for simplicity we shall choose T = 1 and e = k~^, where e is the error 
of the data and k is an integer. The regularized scheme in Theorem 2 can be rewritten as 
the following procedure, where k, ipk, gk are given data and fk is the regularized solution. 

for n from 1 to 20r do 

for j from 1 to 20r do 

1 1 
J j 9k cosh.{zjx) cos{mTy)dxdy 

zj ■= 4r + j; Hj := - ° ° ^ ; 



end do; 

for m from to r do 
coef := 0; 

for j from 1 to 20r do 
w := 1; 
for i from 1 to 20r do 

(TTiTr) z 

if {i 7^ j) then w := w * ^ ^ end if; 

Zj - Zf, 

end do; 

coef := coef + w * Hj; 
end do; 

fk '■= fk + K{m,n) * coef* cos(m7rx) * cos{mry); 
end do; 
end do; 

In the two following examples, we will show the ill-posedness of the problem and how 
our regularization treat it. 
Example 1. Consider the exact data 

ipo{t) = Ti'^e'^'' \ go{x, y) = {1 + cos(7rx)) cos(7r2/). 

Then the system (1) has the exact solution 

Uo{x,y,t) = e"^'^ * (1 + cos(7rx)) cos(7r?/), 
fo{x,y) = —3cos{ny) — 2cos{'Kx) cos{TTy). 

For each integer A; > 1, take the disturbed data 

Vk{t) = ipoit),gk{x, y) = goix, v) + t ^^"^^{knx) cos{k7ry). 
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Then the system (1) has the disturbed solution 



n 



Uk{x,y,t) = UQ{x,y,t) -\- —e '^^ ^ sin^i^kirx) cos{k'Ky), 

K 



TT 



fk{x,y) = fo{x,y) + -{{5k^ -4)8111^ {knx) -2k^)cos{kny). 



We see that 



\9k gWi^n) - ■^' 



fk — /o 



n / 48 

= -j27A;2-56 + -. 
L^n) 4 V k^ 



Hence, if k is large then a small error of the data will cause a large error of the solution. 
Thus the problem is ill-posed and a regularization is necessary. 

Corresponding to the error of the data e = k~^ = 10~^, our regularized procedure 
produces the regularized solution 

fk{x,y) = -2.999721 cos(7r|/) - 1.997145 cos(7ra:) cos(7r2/). 

The error between the regularized solution and the exact solution is 

ll/fe-/ollL^(n) = 0.001441. 

In this case, the exact solution has form of a truncated Fourier series, and hence the 

approximation is very good. 

Example 2. Corresponding to the exact data 

ipo{t) = e\ go{x) = (xcos(l - x) + sin(l - x) - 1) {2y^ - 3y^), 

the exact solution of the system (1) is 

uo{x,y,t) = e* (xcos(l — x) + sin(l — z) — 1) (2?/^ — 3|/^), 
foix,y) = {2xcos{l-x)-l){2y^-3y^) + 

— (2;cos(l — x) -\- sin(l — x) — 1) {12y — 6). 

For each integer k > 1, the disturbed data 

TT 

V^fc(t) = (poit),gkix, y) = goix, v) + t sin^(A;7rx) cos{2Txy) 
corresponds the disturbed solution 

TV 



Uk(x,y,t) = uo(x,y,t) + —e^ sin^ikiix) cos(2TTy), 

k 



TC 



fk{x,y) = /o(x,y) + -((4fcV + 47r2 + l)sin2(A;7rx)-2A:V)cos(27ry). 



We see that 

\\9k — 5'llLi(f^) 

fk — /o ^^ ^ 

L2(n) 



1 

k 



0, 



TT / ,,, , , 487r4 + 247r2 + 3 

-\ 16n^k^ + 327r4 + Sn^ + -— > +cx) 

4 V A;2 
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as A; ^ +00. Thus a small error of the data causes a large error of the solution. 




Q 



Figure 1. The exact solution fo{x,y). 








Figure 2. The disturbed solution /^ with k = 100. 
Now, we use our regularized procedure for e = n~^ = 10^^. The regularized solution is 

fk{x,y) = 0.040435 + 0.426992 cos(7rx) 

-0.431701 cos{ny) - 0.800509 cos(^a;) cos(7r2/). 
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Figure 3. The regularized solution /^ with k = 100. 

The error between the regularized solution and the exact solution is 
II /fc — /o|Il2(q) = 0.059997. To see the effect of the regularization, we note that the 
disturbed solution corresponding to A; = 100 causes a so large error 



fk — /o 



L2(f7) 



1.24 X 10^. Figure 1, Figure 2 and Figure 3 give a visual comparison 



between the exact solution, the disturbed solution and the regularized solution in the 

second example. 
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